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Abstract 

A non-covariant but approximately relativistic two-body wave equation (Breit 
equation) describing the quantum mechanics of two fermions interacting with one 
another through a potential containing scalar, pseudoscalar and vector parts is pre- 
sented. After expressing the sixteen component two-body wavefunction in terms of 
a radial and an angular function by means of the multipole expansion, the initial 
equation can be reduced into a set of sixteen radial equations which, in turn, can 
be classified in accordance to the parity and the state of the wavefunctions involved. 
The adequacy of the reduced equations in describing real problems is discussed by 
applying the theory to QCD problems and the calculation of the energies of bound 
states of quark-antiquark systems is performed to order a^. We show that bound 
states of heavy quarks can be described adequately by the Breit equation for a funnel 
interaction between the particles. 

PACS numbers: 03.65.Ge; 03.65.Pm; 12.39.Jh; 12.39.Ki. 



1 Introduction 



Over the years a number of equations has been introduced in an effort to describe ade- 
quately the relativistic dynamics of a system of two interacting spin-| particles (IQ-IQl)- 
The problem is to determine how the two fermions behave when they are influenced by 
their mutual interaction. The simplest relativistic equation for fermions, the Dirac equa- 
tion, which can describe the quantum mechanics of a single fermion, is not useful for a 
two-fermion system unless the mass of one of the particles is much larger than that of the 
other particle. One equation that has been extensively used, in the past, is the equation 
Breit proposed in 1929, which describes the interaction between two electrons [§, ^ |] 

[E - H^,) - H(2) - Vint{\x^ - ^2|)V^(^1, ^2)] = (1) 

where H(i) = Ei ■ pi + (3m (z=l,2) is identical with the Dirac Hamiltonian of the i particle 
and E is the total energy of the system. The interaction between the two particles equals 
yint = Vc + Vb and it is partly due to an instantenous (static) Coulomb interaction 
and partly due to effects prescribed by quantum electrodynamics described by retardation 
terms 

Vs(r) = --Vc(r)(ai ■0.2 + 0,1- fa2 ■ r) 

which is referred to, as "Breit Interaction". Oi are the Dirac matrices, r{= Xi — X2) is the 
relative position vector and f is the unit vector. 

Although the initial Breit equation involved two electrons 0, we aim to generalise it 
for any two fermions. There are two requirements that the Eq.(|l]) clearly satisfies: (i) 
In the limit of negligent interaction between the particles, Eq.(|l|) implies that the total 
energy of the system equals the sum of the energies of the particles and the stationary 
wavefunction ip{xi, X2) is the product of the wavefunctions ip{x2) which are simply 

the solutions of the Lorentz-invariant Dirac equations for each particle, separately, (ii) In 
the non- relativistic limit, the Breit equation reduces to the Schrodinger equation for a 
two-particle system. 

If retardation effects are ignored, the Breit equation reduces to the Dirac equation in 
the infinite limit of the mass of one of the constituent particles and then, either equation 
can be used. Certainly, the former would provide a more accurate description because 
corrections due to the motion of the heavier particle should be taken into account. 
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Unlike the Dirac equation, the Breit equation lacks a very important ingredient char- 
acterising the relativistic equations, the covariance. Eq.(|l|) is not written in a covariant 
notation since individual terms in the part representing the interaction are not Lorentz 
invariant (the potential is not a relativistically invariant quantity). Besides this, not only 
does the wavefunction in the Breit equation depend on the positions Xi and X2 but also 
depends on one time- variable (rather than an one time- variable each for the two particles), 
a fact that does not allow a Lorentz-invariant formulation. However, it can be consid- 
ered covariant in the centre-of-mass-frame because, in that case, the relative motion of the 
particles is studied and there is only one position vector which comes into the equation, 
namely that one which measures the distance between the two electrons, f. Attempts have 
been made aiming at deriving a covariant Breit equation by means of constraining the 
equations. It is too early to say whether this work is successful 0. 

In spite of the fact that Breit equation manifestly lacks covariance, it appears that 
it can describe efficiently the quantum mechanics of some two-body systems (at least 
approximately to the required accuracy) for long-range interactions with small coupling 
1^, rU|. When the equation was proposed by Breit in 1929 0, QED was known and the 
long range Coulomb type potential Vc{r) was the form of the instantaneous interaction for 
hydrogen-like atoms. In contrast to the static Vc{r) which is the zero-order term in the (v/c) 
expansion of the electromagnetic interaction, the Breit interaction VB{r) is the (v/cf term 
(velocity-dependent) and it constitutes only an approximation to the relativistic interaction 
between the two particles. The study of the fine and hyperfine structure for hydrogen-like 
atoms indicates that the Breit interaction should not be considered on the same footing as 
the Coulomb potential Vj^r). On the contrary, it should be treated as a small perturbation, 
otherwise it does not lead to correct results (|lO|-[0). The fact that the Breit equation 
gives satisfactory results for the Coulomb potential implies that the equation, although not 
covariant, can provide a good description to a two-body system for long-range interactions. 

In this work, we aim to extend the application of the Breit equation to systems of 
fermions that interact with one another through a short-range strong potential. Due to 
the fact that the distance between the two particles is very small, one position vector 
suffices to describe the system which puts in the same footing time and position. Then, 
one might consider the Breit equation is approximately covariant and more compatible 
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with the special theory of relativity, in the case of short-range strong interactions. 

There has been a belief that the internal dynamics of quark-antiquark systems can be 
described satisfactorily by two-body semi-relativistic equations and some attempts have 
been made towards this direction ([]T3|,|T^-[Q). This is a very interesting possibility since 
it suggests a qq bound state could provide a sort of "Hydrogen-atom" for QCD. In this 
paper, we consider fermionia composed from heavy quarks (bottom, charm) and we assume 
that the particles interact with each other through a "funnel" potential |jl6|, Through- 
out this paper, we will handle the Breit equation by taking into account the instantaneous 
interaction between the particles and treating the non-static terms perturbatively. It turns 
out that the Breit equation offers a very good desciption of the systems and this is partly 
due to the fact that the consituents are heavy particles whose speed are very small so as to 
declare that a semi-relativistic treatment gives satisfactory results for certain bound states 
of the two fermions. Unfortunately, this does not seem to hold true in the event of higher 
bound states as well as for bound states of lighter particles because a relativistic treatment 
is necessary. 

One of the advantages of the Breit-equation is that it can be handled relatively easily 
by means of some standard techniques (i.e. multipole technique [^, 0). The 

wavefunction describing the system is a sixteen component function and it can be used to 
obtain radial equations. 

This paper is organised as follows. We start in Section |^ with the introduction of the 
Breit equation and we assume a general interaction between two spin-1/2 particles including 
a vector, scalar and pseudoscalar particle exchange. Using the multipole technique, we 
separate the equation into its radial and angular parts and derive sixteen radial equations 
grouped according to the values of the spin of the system. Then, we present the non- 
relavistic type of the potential which describes the interaction between two quarks [Section 
^ and we derive the radial equations that apply in that case [Section In Section ^, 
we solve the equations for two types of quarkonia, bottomium and charmonium, and we 
calculate the energy levels of various bound states. This consideration allows us to check the 
validity of the Breit equation by comparing the results derived from the theory with well- 
established results for bottomium and charmonium. Section ^ contains brief concluding 
remarks. 
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2 The Breit equation and its reduction to radial equa- 
tions 



We consider a two-body Dirac equation (Breit equation) [0, |, H 

E - 7«(7(^) . p + mi) - 7^ (-7^'^ ■ P + m^) - VU\^)\ ^^(r) = (2) 
describing a system of two spin- 1/2 particles of mases mi and m2, in the centre of mass 



frame, interacting with each other through a static central potential of the form [|10|, p8 

V^nt{r) = Vs{r) + Vp{r) + Vv{r) (3) 
Vs{r) = -(7?^®7f )5(r) 

Vv{r) = +[(7S'^i^^)®(7f7f)],=o^W 
where r = xi — X2, p = Pi = —p2, r = \xi — X2\, and Vs{r), Vp(r) and Vv{r) are the 
parts of the interaction with scalar, pseudoscalar and vector Lorentz structure, respectively. 
The choice of the combination of the 7 matrices which leads to the appropriate Lorentz 
structure is not unique, however, the selection is based on some conditions which should 
be satisfied. To be more specific, 70^'' -|- 70^'' has a scalar Lorentz structure, as well, and 
although it seems reasonable because it couples the potential S{r) directly to the mass of 
each particle (a feature of the scalar potential in contrast to the vector potential which 
couples to the charge of the particles) does not flip the helicities of the fermions in the case 
of the chiral representation. 

The signs in front of the various forms of the potentials can be justified as follows: in 
QED, the vector exchange is actually the Coulomb interaction between particles and if 
their charges have the same sign, then Vy(r) is positive. On the other hand, for scalar 
and pseudoscalar exchange the sign of propagator is opposite to the Coulomb part of the 
photon propagator. 

The superscript i {i = 1,2) which appears in the 7 matrices refers to the particle. The 
reason why in the case of potentials with vector Lorentz structure only the 7g^'*, 70^'' contri- 
butions have been taken into account, is that the 7''*^ matrices give a non-static character 
to the potential since they introduce velocity terms [|. 7^^''(7''^^) is the Dirac matrix 7 
^We recall the relation -y^*-* = ^^^0.^^^ = lo^'^hr ' "^ti^re ■iT^*-' is the velocity of the i particle. 



acting in the subspace of the spinor of particle 1 (2) and it acts on ip from the left (right) 

7(2)^; = ^(7(2))'^ 

The total spin of the system of the two particles is either S = or S* = 1, therefore the 
total angular momentum of the system is either j = I (for S = 0) or j = (for S = 1). 

The form Eq. (^ acquires depends on the representation of the 7 matrices. In the next 
subsection, Eq. (^ is writen in the Dirac-Pauli representation, where 7q*^ matrices are taken 
to be diagonal. This choice will allow us to examine the connection with the non-relativistic 
limit. 

2.1 Dirac representation 

In the Dirac representation, 7q*'' matrix is diagonal 

The spinor ilj{r) is a sixteen-component wave function and it can be represented as a 4 x 4 
matrix 

m = ..,,,„) = ( J- ) (6) 

where the indices +, — are the eigenvalues (+1, —1) of the Dirac matrices 70^^70^"* in the 
so-called double Dirac representation [|^, The left index refers to the first particle 

and the right one to the second particle. By inserting (§) and (^ and (|^) into Eq.(||), the 
Breit equation takes the form 

In order to simplify Eq. (^ and bring it to a form which can be handled easily, we introduce 









^(1) 

cr^ ^ • 






-.(2) 
^(2) 













the following components ||36|, ^| 



I = Po^(V^++T^-) 

I = i(a^«-a^(^))Pi^(^+_±^_+) 
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X 
X 

x' 



)Pi^(^++±^__) 



1/-(1) -►(2) 



(9) 



where 



Pi 



^(1) ^(2) 



(10) 



are the projection operators on states with total spin S* = and S* = 1, respectively, 



Pol state) 
Pi I state) 



I state S = 0) 
|state S = 1) 



(11) 



The components 0, , x, X and </), 0°, x, x° correspond to spin S* = and 5 = 1, respec- 
tively and they are functions of r. Noting that 



and by means of the identities 

p |'^(1)^{2) ^(1)^(2)\ 



J2)^ 



—2SikPo 
iPoeikii^^'^ - 



2Po5. 



ik 



(12) 



(13) 



Eq. (^ leads to the following set of component wave equations |jTO|, ^ 



E + S{r) 


- P{r) 


- V{r) 


<f- 




= 


E + S{r) 


+ P(r) 


-V{r) 


0- 


(mi+m2) ^0 


= 


E - S{r) 


+ P(r) 


- V{r) 


x'- 




= 


E - S{r) 


- P(r) 


- V{r) 


X- 


(mi-m2) 
2 


= 


E + Sir) 


- P(r) 


- V{r) 


X- 




= 


E + S{r) 


+ P(r) 


-V{r) 






= 


E - S{r) 


+ P(r) 


- V{r) 


0- 




= 


E - S{r) 


- P(r) 


- V{r) 






= 



(14) 
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The type of the potential we consider is central therefore the next step we will follow is 
the introduction of a method that will eventually separate the angular from the radial 
dependences as it was performed in the three dimensional Schrodinger equation. To this 
end, it is convenient to introduce the derivative 

^ {Sik-niUk)^ (15) 



where n = is a three dimensional vector, while lies in the perpendicular plane 
UiQ^ = ni{5ik — nink)-^ = ^k-^ — "^k^^ = . Some useful identities that are going 
to be used are the following 

dr rsn^ dr 

-4--n = {5ik-nink)^ni = {5ik-nink)5ki = ^-l = '2 



Then, = ^ = (-^ + n-^f = ^ + + This form of the Laplacian 

indicates that -^=pr was correctly regarded as the angular part of the Laplacian. 

2.2 Expansion of the wavefunctions 

The component functions defined by (|^) and (§) are classified in two groups: the scalar 
(0, X°) which refer to states with 5 = 0, and the vector (</>, </>°, X, X°)) which cor- 



respond to states with 5 = 1 |jTO|, ^ ^ These functions can be expanded in the 
following way: 

(i) . Scalar functions (S* = 0). 

The scalar functions can be written as 

B{f)=Y.B{r)Y,^{ri) (17) 

where B{r) stands for any of the four 5 = component functions, B{r) is the radial part 
of the function and Y^„i(n) are the spherical harmonics depending only on the angles. 

(ii) . Vector functions (5 = 1). 

We define an operator Sk which acts on any of the four vector components [(H), (H)] as 
follows 
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— » — * 

where (...) equals + ip-+) (for 0), ('?/'+- — 4'^+) (for 0° ), ('?/'++ + ) (for x); 

— * — * 

('?/'++ — ) (for x^) s-iid A(r) is any of the four vector components. We notice that 

— * 

the operator 5*^ does not act on A(r) from the left since in that case it would give zero 
(^i(CT^^) + o'(^))4(ct(^) - <T^^^)iPi §^ 0^. From the definition of Sk{A{r))\i, it is clear that 
Sk{A{r))\i = ieiki (pi. By making use of the multipole technique |jlO|, ^ it is possible 

— * 

to expand any vector component A{r) into three parts: "electric" (Ae), "longitudinal" 
(Ai) and "magnetic" (y4m)defined as 



Amir) = ^Amir)Yj^{n) 



r -> " V ' " V ' 

Ai{r) A^{r) 

(In the rest of our study, for the sake of simplicity, we will the drop the summation symbol 
y^). It may seem that j in Yj^ is the orbital momentum rather than the total one. But this 

j,m 

is not true for the function A{r) because the vectors to which it is proportional, depend 
on angles. Actually, if S and L are the total spin and total orbital angular momentum 
operators, respectively, then for each part of A{r), for instance Ae{r) = nAe{r)Yjm{'n) , 
one has 



J,(A^(r))|, = S,{A^{r))U + L,{A^{r})U = 

d 

= ieikiniAe{r)Yjm{r) - iekiani — {niAe{r)Yjm{n)) 
dY-min) 

= -iekianiUiAeir) — = niLk{Ae{r)Yjm{n)) = 



dua 
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J A,ir} = nL\A,{r)Y,mifi)) = j(j + l)Ae(fO 



which implies that j is the total angular momentum. By inserting [([T7| ) ,([T9| ) j into Eqs.([T^), 
we manage to eliminate the angular dependences and the Breit equation reduces to the 
following set of sixteen radial equations 
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E + S{r) 


— P(r) 


- V{r) 




- 2 0-(^ + ^)0e - 


1 , 

70« 


= 


E + S{r) 


+ P(r) 


- V{r) 




f?77,i -l-?T),0 '1 1 n 

2 ^ 




= 


E - S{r) 


+ P(r) 


- V{r) 





("tTIi — TTlo'l / f] , 9\ 

2 X (,; + f)Xe 


I 


= 


E - S{r) 


- P(r) 


- V{r) 


X- 


2 X 




= 


E + S{r) 


— P[r) 


- Vyr) 




2 Ae "T 




= 


E + S{r) 


- P{r) 


- V{r) 


Xi - 


2 A/ A 




= 


E + S{r) 


— P(r) 


- Vyr) 




fmi +m9) n 
2 




= 


E + S{r) 


+ P[r) 


— V[r) 


n 


fmi -{-mo) 1 1 if) 
^ 2 Xe+,0m 




= 


E + S{r) 


+ P[r) 


— V[r) 


n 

X/ 


2 Xl idr + r)0m 




= 


E + S{r) 


+ P(r) 


- V{r) 


n 

Am 


( Tn,-i ~\-Tn,o^ . 7^74-1^ /fi 
2 Am 1 r 


ii + ^)0° 


= 


E - S{r) 


+ P(r) 


- V[t) 


0e - 


(mi— m2) iO 1 q! J,0 
2 dr^ 




= 


E - S{r) 


+ P(r) 


- V{r) 


0^- 


{mi-m2)^0 i(i+l)^0 




= 


E - S{r) 


+ P(r) 


- V{r) 




(mi-m2) jO 
2 r'm 




= 


E - S{r) 


- P(r) 


-V{r) 




- ^^0e + ix^ 




= 


E - S{r) 


- P(r) 


-V{r) 


0? 


- ^^^0^ - (f + ^)X^ 




= 


E - S{r) 


- P(r) 


-V{r) 


0^1 


(mi-ma) , , j(j+l) i 
2 Vm 1 J, Ae ~r 




= 



(20) 



At this point, we will concentrate on the properties of the scalar and vector components 
as well as the new components to which they have been expanded. As pointed out earlier, 
the scalar components 0, 0*^, x, describe the S = states, while the rest, the vector 
components, describe the states with = 1. In the latter case, for x, there are three 
states characterised by the spectroscopic signatures ^(j — 1)^,^ (j + l)j and ^jj (in the 
atomic notation '^^'^^Lj) described by the wave-functions xi=j~i,Xi=j+i Xi satisfying 



the following relations Bq, 37 



Xe 
Xl 
Xm 



2j+l 



Xi-- 



-1 + 



2j+lXl=3+l 



V j(j + i) (-v^Xf=j-i + 4^iXi=]+i) 

^/j(J+l)Xi=j 



(21) 
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Certainly, this can be generalised for the rest of the S = 1 components: x^, (f) and 
(p . The "magnetic" component is chosen to have I = j while the "electric" and the 
"longitudinal" components have / = j ± 1 mixed if j > (and / = 1, if j = 0). If 
j = 0, then there should be only one x component and we note this is satisfied since 
Xi = Xm = and Xi=j-i = Xi=j = 0, Xi=j+i = Xe- There are only two states with j = 0, 
'So (j = 0,/ = 0,5 = 0) and ^Po {j=OJ = l,S = 1). 

The expressions (H) and (||) indicate that, for S = 0, the "large-large" {4'++) and "small- 
small" {ip ) components are contained in the x^'X^ and the "large-small" ('?/'+-) and 

"small-large" {ip h) components appear in the 4>,4> ■ In addition, all "electric" and "lon- 

gitudinal" components of the X^X ^4*^4* correspond to states with magnetic quantum 
number = ±1 while the "magnetic" components describe states with = 0. On the 
other hand, the states with 5 = have the spectroscopic notation 'jj and the "large-large" 

and "small-small" {ip ) components are contained in the 0, (f)^ and the "large-small" 

and "small-large" components appear in the X^- AH 0, (f)^, Xi X^ correspond 

to states with = 0. 

The parity of the system equals P = 1)', where 77 = 1 or —1 for fermion-fermion 
or fermion-antifermion system, respectively. The sixteen components have the following 
parity: 

(i). 0, 0° : P = ri{—iy because they both describe the S = states, therefore / = j, 



n 



^e,0°,0h0?: 



P = —1]{—1)^^' (or —r]{—iy ' which gives the same result) = //(—I)-'. 



According to (0) the "electric" and "longitudinal" components are combinations of states 
with I = j + 1 and I = j — 1 and this justifies the exponent j + 1 (or j — 1). The (— ) sign 
in front of rj is due to the fact that all these functions are combinations of "small" and 
"large" components, 

(iii) . Xm,Xm '■ P — v{~^y ((PD implies that the "magnetic" component corresponds to 
I = j states), 

(iv) . X, x° : P = — ?7(— !)•'. Both X) describe S = states therefore / = j, however, in 
contrast to the parity of the 0, 0°, there is a minus sign which is attributed to the fact that 
X, X^ are combinations of "small" and "large" components, 

(v) - Xe, Xe^ Xu X^i - P = -^(-1)^ (as in (ii)), 

(vi) . 0^,0^:P = -77(-l)^ (asin (iii)). 
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Again, the combination of "small" and "large" components accounts for the minus sign 
in (v), (vi). To summarise, the components 0, 0°, 0^ 0?' Xm, Xm have parity P = 

ri{—iy while the components x") Xe, X% Xu X?; 0m! 0m have parity P = — !)•' . We call 
the former case, Pseudoscalar Particle Trajectory (PPT) while the latter is called Vector 
Particle Trajectory (VPT) PB|, |57[. The "large-large" and "small-small" components in 



the PPT are contained in the 0, , Xmi Xm 

and have spectroscopic signature ^jj or ^jj. In 
the VPT, the "large-large" and "small-small" components are contained in Xe,Xe;Xi)X? 
and have signature ^(j — 1)^ or ^(j + 1)^. 

Returning to Eqs. (pOD and by making use of the previous discussion, that set of equa- 
tions can be split into two sets of equations according to whether they belong to the PPT 
or VPT regime. To be more specific, if the two fermions do not possess the same mass, 
the following two sets of equations are obtained from Eqs.(pOD [|1( 

(i). PPT, ij, or3j„P = r/(-iy-. 



36, 37 



1 

2 


E + Sir) 


- P(r) 


-Vir) 


<f- 


-^^^0-(i: + f) 


- ^0« 


= 


1 
2 


E + Sir) 


+ P(r) 


- Vir) 


0- 


(»ni+m2) ^0 




= 


1 
2 


E + Sir) 


- P(r) 


- Vir) 


Xm 


(mi+m2) 
2 Xm 




= 


1 
2 


E + Sir) 


+ P(r) 


- Vir) 



Am 


(mi+m2) 1 i(i+l) 
2 Am 1 J, V 




0? = 


1 
2 


E - Sir) 


+ P(r) 


-Vir) 


0e - 


(mi-m2)j0 1 d iQ 
2 '^e dr^ 




= 


1 
2 


E - Sir) 


+ P(r) 


-Vir) 


0^- 


(mi-m2)^0 i(i+l)^0 




= 


1 
2 


E - Sir) 


+ P(r) 


-Vir) 


0m 


(mi-m2) iO 
2 ^m 




= 


1 

2 


E - Sir) 


- P(r) 


- Vir) 


0°- 


- (^^0e + ix^L 




= 


1 

2 


E - Sir) 


- P(r) 


- Vir) 


0?- 





Am 


= 



m) 
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(ii). VPT. ^^a±i),.p = 


-ni-i 








1 

2 


h — b{r) + r[r) 


- V{r) 


n 


ivr\A — Tn.o\ / /-/ , 9 \ 1 
2 ^ ^dr 7- / Ae ^ At 


= 


I 

2 


E - Sir) - Pir) 


- V{r) 


X- 


2 X 


= 


1 
2 


E + Sir) - Pir) 


- V{r) 


Xe 


2 Xe "T drX 


= 


1 

2 


E + S{r) - P{r) 


-V(r) 


Xi - 


(mi+m2) j(i+l)^.0 
2 X/ r X 


= 


1 

2 


E + Sir) + Pir) 


- V{r) 


xl 


(nii+m2) 1 1 iO 
2 Ae "T ^V^m 


= 


1 
2 


E + S{r) + P{r) 


- V{r) 




- ^^x/ - (f + 


= 


1 

2 


E - S{r) - P{r) 


-V(r) 




- + + (1: + i)x? 


= 



(23) 



On the other hand, in the case of equal masses mi = m2 = m, the set PPT sphts into 
two subsets, one with spectroscopic signature ^jj and a second with ^jj while the VPT 
remains unchanged 

(i). PPT, 'jj,p^n{-iy. 



V{r) 


/-m0-(|: + f)0e-^0/ 


= 


V(r) 


(f) — m(f)^ 


= 


V{r) 


'i>e + 


= 


V{r) 




= 



(24) 



(ii). PPT, ^j^,P = r,{-iy . 



E + S{r)- P{r)-V{r) 
E + S{r) + P{r)-V{r) 
E - S(r) - P(r) - V(r) 
E - S(r) - P(r) - V{r) 



Xm '^Xm 
Xm - 'n^Xm + 

€ + h'm 










(25) 
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(iii). VFT,'ij±l)j,P=-r)i-iy . 



1 

2 


E - S{r) 


+ P 


y) 


- V{r) 


x' 


( (jj. ~l~ ^ ) Xe J. Xl 


= 


1 
2 


E + ^(r) 


-P 




- V{r) 


Xe 


- nixl + 


= 


1 
2 


E + S{r) 


-P 


t) 


- V{r) 


Xl 


- mx? - 


= 


1 
2 


E + S{r) 


+ P{ 


r) 


- V{r) 


xl 




= 


1 
2 


E + S{r) 


+ P{ 


r) 


- V{r) 


X'l 




= 


1 
2 


E - S{r) 


-P 




- V{r) 




+ + (f + 


= 



(26) 



The remaining two components x, 0^ vanish in the event of the particles having the same 
mass m. 



3 A Potential for quarkonium 

In the Introduction, we discussed the lack of covariance characterising the Breit equation. 
The wavefunction of the system of the two fermions depends on the position of the two 
particles (that is, two vectors, each for the two particles, are necessary) while it depends 
only on one time variable. This fact does not allow a Lorentz-invariant formulation. How- 
ever, if a very strong, short-range potential is considered to govern the interaction of the 
two particles, we can assume the two particles are so close to one another that just one 
position vector is enough to describe the motion and the quantum mechanics of the system. 
In this case, the number of the position vector variables equals that of the time variables 
and we can assume that the equation becomes approximately covariant. The question 
rising is whether there exist such interactions between particles and if there are systems of 
fermions to which the Breit equation can be applied and subsequently be tested. Talking 
about small distances, we are led to the plausible question of what happens at distances 
up to 1 fermi which is a feature of the strong interactions between quarks. It is known that 
there are systems, quarkonia, which are bound states of quarks and antiquarks. Actually, 
experimental data about many of these states exist |^T|, therefore, one could maintain that 



they constitute a laboratory where the validity of any theory aiming at describing them 
can be tested. 
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In the non-relativistic limit, the quark-antiquark interaction is dominated by vector and 



scalar potential [|1^, |3^. The former, which is of the Coulomb type consitutes the interac- 
tion at small distances and corresponds to one gluon exchange. In that case, quarkonium 
can be viewed as the "hydrogen atom" of strong interaction Physics. By contrast, at large 
distances, confinement dominates and the interaction is proportional to the interparticle 
distance. The confining part of the potential has scalar structure [jlO|, ^ and thereby the 
interquark potential consists of a Coulomb part which is of vector type and a linear part 
which is of scalar type 

VUr) = Wector(r) + Vscaiar(r) = -^^^[(7?S?)) ® (7^7^)] + ^r^'^ ® 7^ ) (27) 
where is the coupling constant of the strong interaction, is the relevant momentum 



transfer and k is the string constant ||3^ . 

Certainly, the non-relativistic part of the potential is not sufficient if we wish to test 
the Breit equation. The contribution of relativistic corrections are required and we must 
take them into account. At this point, we recall that the original Breit equation contained 
terms called Breit interaction that account for relativistic corrections. The inclusion of 
the Breit interaction in the original equation does not lead to correct results ([|T^-||T^) for 
certain quantum-mechanical systems. Unless this term is treated as a small perurbation, 
the results differ from the well-established data. 

iTTT^ rrr^ 



By studying the scattering amplitude for a particle-antiparticle collision [1^, ^ 



can derive an effective potential in the Pauli approximation of the Breit equation. That 
potential will account for the interaction and it will include both relativistic and non- 
relativistic contributions. We need bear in mind that the Breit interaction is the second- 
order term in the (f/c) expansion of a Coulomb interaction. In the Pauli approximation, 
the vector part of the quark-antiquark potential contributes the following term to the 
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interaction 

t^Vector (P, r) 
V2 



{-|t} + ^l + ^2 + V^3 



2a, 



2as 



3mim2r 

L S 



-a 
P 



r p 



(28) 



3^ 



rS 



'1 



+ 



9mim2 



5ir) 



where the {...} terms correspond to the static part of the interaction. The rest of the 
terms are responsible for retardation effects and they can be associated to the Breit in- 
teraction. We emphasise that this derivation is in analogy with the procedure followed 
towards obtaining the retarded part of the potential in a particle-antiparticle electromag- 
netic interaction (scattering). We recall, that in that case, the contribution of the Di^ 
{i, k = 1,2, 3) part of the photon propagator in the Coulomb gauge was taken into account 

On the other hand, the scalar confining potential contributes only to the static part of 
the interaction. Unlike the vector part of the interaction, the scalar term of the potential 
cannot lead to the emergence of retarded effects since there is not an equivalent of the 
photon propagator that would result in such a contribution. The scattering amplitude 



(29) 



leads to the following interaction operator [Ty, |33 

t^Scalar (^ = {t^r} + V4 

where {...} merely imply that the interaction is characterised by static behaviour. For 
particles with the same mass, V2, V4 correspond to the spin-orbit interaction. Although 
both potentials, (Vvector) and (Vscaiar) are attractive, they give the opposite sign of spin- 
orbit interaction, a feature which is going to result in a quite interesting effect regarding 
the ordering of the states (Section |5]5). In addition, we should stress the absence of a 
spin-spin interaction in the confining potential meaning that tensor forces are contained 
entirely in the Coulomb-like vector potential. The absence of spin-spin and tensor terms 
in the scalar interaction as well as the different form of spin-orbit terms between the linear 
and Coulomb-like potentials constitute the most striking difference between the two types 
of interaction. 
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So far, the form of the potential describing strong interactions has been discussed to 
some extent but nothing has been mentioned about the efficiency of such a potential. The 
charmonium and bottomium spectra are well-known therefore the accuracy of the 
potential besides the adequacy of the Breit equation in describing strong interactions can 
be tested. Next, we will follow a procedure towards calculating the energy levels of various 
bound states of quark- ant iquark systems, however, we have to obtain first the equations. 
In the following section, we will derive the radial equations for those systems as they are 
produced by means of the Eqs.[ (|2^) , (|25|), (|26|) ]. 



4 Application of the Breit equation to quarkonium 



The quark- ant iquark bound states states can be classified in three categories according to 
their spectroscopic signature ||10| : 

(i). States with j = l,S = 0. After making the substitutions V{r) = — S{r) = —nr 
and eliminating all components but (p^, Eqs.(p^ lead to 



+ 



K — 



3 r2 



dr \r 

jU + 1 



— m 



3 r . 



+ 



E + - — I — K r 



3 r 



/(r) = 



(30) 



(n). States with j = I, S = 1. 



Similarly, Eqs.(l25|) give 



dr"^ dr \r E + nr + "^"^ 



4 Os 
3 



+ 



m 



3 r / I 

^ 3 



i5 + l^V_^V2 

3 r 



J(J + 1) 



3 r 3 r 



X^(r)=0 



(31) 



(ni). States with j = / ± 1, S* = 1. 

Finally, Eqs.(EB) can be reduced to the following two coupled differential equations 



dr'^ 



+Xe(r) 



dr 

{e + 



3r2 ^ 



2 

r E + ^ — K,r 

3r 



4«3 



+ 



-2k (e + 



3r 



Kr] 



+ 



4aa 
3r 



+ Kr 



4 V 3r 



Kr . 



m 



E + 



4qs 
3r 



Kr 
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= X%r) 



-2k {e+^ 



ii 



+ 



E- 



4as 
3r 



— ACr ) 1 



(32) 



X 



1 
4 



3r 



4a, 



3r2 ^ 



2 

2 



X 



K,r 



m 



j{j + 1) 



1 -^ + ^ 1 

3^:2- + « 



+ 



K,r 



(33) 



The behaviour of the component wavefunctions in the above equations near the origin can 
be deduced easily if we notice that at very small interparticle distances, the dominant part 
of the potential has a vector, Coulomb-like form. Thus, the wavefunctions for small r look 
like 



state 








, 7 = 


state 


% , 






, 7 = 


state 


3/ 


f rxl 
\ rx^i 




, 7 = 



4 9 

-cy 

9 s 



-1 + JjU + 1) 



4 2 



j(j + l) + l-|«2 



(34) 



-1 + Jj(j + 1) 



4 9 

-Q; 

9 s 



On the other hand, due to the inclusion of the confining potential, the wavefunctions have 

1 2 

a different behaviour at large distances. They look like ~ e~*'^'^ . The behaviour of the 
wavefunctions near the origin and at large distances will be used in the next sections in 
order to solve the differential equations describing the dynamics of the quantum mechanical 
systems. 



5 Solution of the equations 

As pointed out in the previous sections, the energy levels of the various states are derived 
from the solution of the above radial equations, however, those values do not represent the 

complete energy since the potential used is merely the static potential. The contribution of 
the retarded part of the interaction should be included to obtain the relativistic corrections 
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to the energy. The static part consists of the {...} terms in the expressions [(|28|), p9|)] 
while the rest are the retarded terms. At this point, we should recall that we do not expect 
satisfactory results unless the Breit (retarded) terms are treated by first order perturbation 
theory for particular states. The application to QED problems (|]TU|-[]T3[) demonstrated 



this assertion is correct. In that case, the total energy of the system is given by the 
expression 

-Etotal = -^static + (state | Retarded | Statc) (35) 

where |state) stands for the stationary unperturbed states for the Coulomb potential. The 
same procedure can be followed when the potential has the form (^7D, however, unlike 
QED problems, the eigenstates are not known and they should be calculated before we 
proceed. Unfortunately, an exact solution for such a potential cannot be obtained and 
an approximate method should be tried. One way is by treating the Coulomb term as 
a small perturbation [O. Another method, which will be employed in this paper, is 



to use the three-dimensional isotropic harmonic oscillator (TDIHO) eigenstates |T0|, P3[ . 
The potential of the harmonic oscillator provides confinement of quarks and the important 
feature of this potential is that the wavefunctions and the matrix elements can be calculated 
easily and they can possess an explicit form. The parameters of the wavefunctions will be 
adjusted so as to fit the numerically calculated eigenfunctions for the potential (^Tf) . 
The interaction potential of TDIHO equals 

Vtdiho = ^Atw^'"^ (36) 

where /i and uj are fitting parameters. The energy levels of the system are given by the 
expression 

E{nr,l) = (2nr + l + ^u, n^, / = 0, 1, 2, ... (37) 
and the normalised radial parts of the wavefunctions which are going to be used are the 
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following 

Ris = 2 (^) exp (-^A , R,P = ^rB,p{r) 



Bisir) 

R2S = VQ{l-lX2sr^)B2s{r) , R2P = {5 - 2X2pr^) B2pir) 

Ess = ^i^5- 20X,sr^ + AXlsv') B.sir) , i?4S = 7= (105 - 210X,sr'+ 

+S4Xlsr^-8Xlsr^)B^s{r) 
= (945 - 2520A55r2 + l512Xlsr'- R,s = ^gg^eoo (^0395 - 34650X,sr'+ 

-288Xlsr^ + IGXtsr^) B^sir) , +27720Xlsr^ - 7920Xlsr^+ 

+880Xlsr' - 32A^5riO) B^sir) 

where -Bindex = B {Xis Aindcx)- Xis = X = fiu while the rest of the As are functions 
of Ai5 and their values will be determined from experimental data. This is done because 
the harmonic oscillator is an approximation of the "funnel" potential [|l^, ^ therefore 
the values of the various parameters need to be adjusted to be in agreement with the 
experimental values of the corresponding quantities. It is necessary to emphasise that the 
As have a mass dependence implying that the relation among them which will emerge is 
not the same for all quarkonium systems. The parameters /i, u, in turn, can be estimated 
by comparing the values of the quantities in which they enter with those derived from 
experimental data. The leptonic widths, the mass differences and the various energy levels 
of some certain states are quantities that can help us to estimate not only those parameters 
but also other parameters such as the coupling constant a^, the string constant k, and the 
quark masses. To achieve this goal, we will take advantage of the knowledge of the energy 
levels of the charmonium and bottomium systems |^T]] . 



5.1 The Bottomium system 

We aim to estimate the energies of the twelve states of the bottomium system which are 
spin triplet states, i.e. S = 1 {ortho-bottomium) . The 1 ,0+^,2++ states are vector 
mesons and the static part of the interaction between the constituent particles leads to the 
description of their dynamics by Eqs.[(|32D, (pS])]. On the other hand, the rest of the states, 
denoted with are pseudoscalar mesons and they are described by Eq.(^). Finally, six 
more pseudoscalar meson states (Table ^ (with J^*" = ^, 1+^) will be considered and 
Eq.(pnD describes those states. Unlike QED, the explicit form of the potential of QCD and 
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the coupling strength are not known therefore a method towards estimating them should 
be formulated. The procedure becomes more difficult if we remember that the solution of 
the above equations is not enough to get the energy spectrum since the values calculated 
account only for a part of the whole energy. The rest of the contribution results from 
relativistic effects due to the non- instantaneous interaction which should not be ignored. 
The combination leads to the problem resisting an analytic treatment and only a numerical 
method seems capable of solving it. 

Since the coupling constant of the system as well as the mass of the bottom quark 
rrih are not known, the use of some of the states to determine these values is inevitable. In 
addition, there are some more parameters which need to be estimated: the string constant 
K, and the As. Before we embark on calculating the various quantities, we will first attempt 
to express the energies of the states in terms of the unknown parameters. 

The energy of the system consists of two parts, one coming from the non-perturbative 
solution of the Eqs.[(pOD-(p^)] and another derived from the perturbative treatment of the 
non-retarded parts of Vi, V2, V3 in ( [28| ) by repacing mi, 1712 with rrib = m. The former will 
be called -Estatic while the latter which is non-instaneous, in nature, will be called i^non- static- 
The combined result in the total energy of the system -Egystem is 

/ 2a, 



system 



+ {Vi + V2 + V3) = E, 



static 



-^non — static 



3m2 

, 2 



2 cP 2d 1-2' 

r dr'^ dr 



+ 



9m^ 



(38) 



S12 



The spin-orbit perturbative correction and S12 vanish for ^-states while the 5 functions 
give a non- zero contribution. By contrast, the P-states are characterised by the opposite 
behaviour (l ■ 5, ^ 0, {5 (r)) = 0) . In Table ^ the energies of various states are sum- 
marised and we notice that the knowledge of the As is necessary to obtain the values of 
the corrections due to the Breit terms. The decay widths of the states will be used, next, 
to evaluate the relation of the As and other parameters to some measured quantities. The 



"Corrected Van Royen-Weisskopf" formula [[18 



V(V 



e e 



1^.(0)1^ (1 - ^a.) 



a 



137' 



^-3 



(39) 
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relates the leptonic width, i.e. e"'"e~ decay of the neutral V vector mesons, and the wave- 
function of bottomium at the origin, rriv stands for the mass of the meson. Although the 
coupling constant has not been estimated yet, the fact that we are dealing with strong 
interactions suggests that the radiative corrections are so large that only ratios such as 



r{V'/V) = 



T{V' 



e+e or 



|^'(0)f 



e+e or yU+yU 



(40) 



can be calculated reliably because the corrections are suppressed and do not appear in the 
above expressions. In (^OD, V' is another vector state of the quarkonium system having 



mass m'y. All 1 vector states can decay into an electron-positron (or muon-antimuon) 
pair and the five ratios of these decays with respect the "reference decay" of the T(IS') 
state are given from experiment |1T[] 



r{2S/lS),, 

r{AS/lS) 
r{5S/lS) 
r{6S/lS) 



rp2 


2 

V'T(25)(0) 


rp2 

-^T(IS) 


2 

V'T(1S)(0) 

2 

^T(3S)(0) 


-^T(3S') 
rp2 

-^T(IS) 


2 

^T(1S)(0) 

2 

V'T(4S)(0) 


-^T(4S') 
rp2 

-^T(IS) 


2 

^T(1S)(0) 

2 

^T(5S)(0) 


rp2 

-^T(55) 
P2 

-'^T(IS) 


2 

^T(1S)(0) 

2 

^T(6S)(0) 


F2 


2 

^T(1S)(0) 



0.4 ±0.1 



0.37 ±0.06 



0.18 ±0.04 



0.23 ±0.05 



0.10 ±0.02 



(41) 
(42) 
(43) 
(44) 
(45) 



(the index /i in (^Tf ) and (^) indicates that the vector mesons involved decay into muon- 
antimuon pairs). Each one of r{nS/lS) {n = 2, ...6) is a function of the energy of the 
system and two As, i.e., Ai^ and Xns- 

The relations [(^lD-(^5D] can be used to establish the relation among the As, energies 
and m. (|38| ) is another useful expression which will be employed in order to calculate the 
energies of the various states as well as the values of a^, k, m. The part of the energy due to 
the instantaneous interaction should be combined with the non-retarded corrections so that 
the energy obtained is in agreement with experiment. This is called "fitting procedure" 



21 



and it will help us to determine the four unknown parameters as, K,,m, Xis- The rest 
of the As, associated with the ^-states, are related to them through the decay widths 
r {V ^ e^e~) whose values are well established. There are two more As, namely, Aip 
and A2P which cannot be calculated accurately by means of the transitions decays, as 
pointed out in the previous paragraph. The total number of the unknown parameters, 
i.e. six, requires six energy values to be set as input values. We choose the energies of 
the T(l^), T(2^), T(3^), T(4^), Xw(l^), Xbi{2P) states as the input values. Although P- 
states with the same principal quantum number can be used, we avoid it because they 
correspond to energies very close to one another which may lead to wrong results. The 
harmonic oscillator wavefunctions constitute an approximation to the real wavefunctions of 
the funnel potential and in the stage of optimising them, the method may be too sensitive 
to small energy differences. Thus, we choose the following values established by experiment 



El 

E[T{1S)] = 9460.37 ± 0.21 MeV 

E[T(25)] = 10023.30 ±0.31 MeV 

E [T{3S)] = 10355.3 ± 0.5 MeV . . 

£;[T(45)] = 10580.0 ± 3.5 MeV ^ ^ 

E [xbi{lP)] = 9891.9 ±0.7 MeV 

E[x6i(2P)] = 10255.2 ± 0.5 MeV 

The energies of the T states will help us to determine the relations between the Xns 
{n = 2,3,4) and Xis through the ratios [(|?T|)-(|D]. 

The mass of the bottom quark is regarded as a parameter which will be evaluated, 
however, its value will be contstrained in the range 4.1 GeV < nib = m < 5 GeV (running 



mass) [^]. The energies due to the static interaction can be expressed easily in terms 
of the total energy and the parameters. From (|3|), -Estatic = -Etotai - -^perturbation which 
means that the numerical solution of Eqs.[(pT|)-(|33|)] can lead to the determination of the 
parameters, provided that the relations connecting them dictated by the corresponding 
quantities are satisfied. Thus, the main problem reduces to the solution of the differential 
equations. Certainly, it is important to observe that Eqs. [([30|)-(|33|)] exhibit a singularity 
at 

IE [T7e\^ 4^ 
This is a feature which is not present in the initial Breit equation but it emerges after 
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the reduction of it to radial equations. Although the equations are singular at r = as 
well, due to the centrifugal term ''^''^^^ {j is the total angular momentum of the system), 
the new singularity implies the appearance of a turning point that is energy dependent 
Pm , |T3j, The question rising is whether the emergence of this singularity causes any 
problem to the solution of the differential equations. Due to the fact that we consider only 
a short-range potential and the wavefunction tends to vanish at distances (r > 1 fm), we 
will not experience any difficulty in solving the differential equations. It will be shown that 
tq is bigger than 10 fm (for the lowest energy) which implies that the singularity does not 
cause any problem in the range where QCD applies, on the grounds that the equations 
have meaning in this particular range. 

The X cind T states of the "fitting procedure" are described by Eq.(|31|) and Eqs.((|3^), 
(0)), respectively, and these equations can be reduced to a set of twenty, first order dif- 
ferential equations. The problem is a boundary-value problem and two boundary points, 
-Rinitiai, -Rfinai uccd to be Specified between which the integration will take place. We choose 
i?finai — 3 fm where the wavefunction and its derivative are taken to vanish due to the 
confinement. On the other hand, the choice of -Rinitiai requires more care. At small dis- 
tances, the Coulomb potential dominates the interaction and since the potential is singular 
at r = 0, we are not allowed to start the integration at that point, therefore we choose 
a point which is much smaller than = ^ (A: QCD scale parameter). The QCD scale 
parameter A is not an independent parameter, it is related, instead, to k Q through 
the expression 

SttA^ 

K = 48 

27 ^ ' 

We set -Rinitiai = 10^^°rc. The wavefunction behaves as r''' at r = -Rinitiai while it vanishes 
at r = -Rfinai- According to (|34D, 



^ ^ -1 + VjU + 1) - |«s , for Xbi states ^^^^ 

+ 1) + 1 - ftts , for '^^ Xbo, Xb2 states 

The reason why the bigger value of 7 was preferred in the ^li±i states is because the 
wavefunction goes to zero faster. 

The numerical integration of the differential equations is achieved by means of a Runge- 
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Kutta-based routine and it leads to the following values for the parameters 



ITOl 



as = 0.36 ±0.04 

K = 0.793 ±0.008 ^ 

tin 

mb = m = 4.987 ± 0.009 GeV 



(50) 



Xis ~ 0.829 (GeV)2 

Aip ~ 0.107 (GeV)2 
Asp ~ 0.120 (GeV)^ 



which, in turn, give 




(41) 



0.147 (GeV)2 
0.126 (GeV)2 
0.072 (GeV)2 
409 ± 40 MeV 



(0) 
(§) 



(51) 



Now that the parameters have been estimated, they can be used to calculate the rest 
of the energy levels of the system. The states that are going to be considered can be 
classified in three groups according to the set of equations that describe their instantaneous 
behaviour or, equivalently, the spectroscopic signature: (i). ri{lS),ri{2S),ri{3S),ri{4:S), 



h{lP),h{2P), with k, (ii).X6i(lP),Xw(2P), with =^ lu and, finally, 

(iii). T(15), T(25), T(35), T(45), T(55), T(65), Xm{lP), Xm{lP), Xmi2P), X62(2P), with 



"^^^^Lj =^ li±i. By applying Eqs.[(^0|)-(^)] to the cases (i), (ii) and (iii), and by imposing 
the same boundary conditions, as previously, the energies of the states can be obtained. 
The As of the wavefunctions entering the equations are dependent on the energies of the 
corresponding states, therefore the expressions [(0), (|45|) ] should be employed. In addition. 
Table |l] and Table provide the contribution of the Breit terms to the energy of the 
system. If all these elements are taken into account, the energy values of the above states 
are calculated and they are summarised in Table ^ and Table ^. The experimentally 
estabhshed energies are also stated. 

The two remaining As, namely, As^, Xqs are calculated and they are equal to 



(i). we notice that the states predicted from the theory have the right accession and almost 




0.082 (GeV)2 
0.045 (GeV)2 



(52) 



Looking at Table ^ and Table ^, there are some points which need to be emphasised: 
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the correct energy difference between the P states emerges. Actually, there is not a very big 
discrepancy between the theoretical results and the experimental data as it can be deduced 
from Table |^. For higher states, it is getting bigger and the discrepancy becomes larger 
and this is quite reasonable since those states are characterised by large speed, therefore 
the f ^ approximation to the Breit terms constitutes a rather rough approximation, 
(ii). the part of the energy due to the instantaneous interaction accounts for more of the 
energy, however, the contribution of terms due to retardation is very important to acquire 
the correct splitting of the P states. At this point, an interesting feature of the funnel 
potential should be mentioned. The spin-orbit interaction coming from the Coulomb-like 



part of the potential (28) differ from the spin-orbit interaction resulting from the linear 

part ( p9D in the sign. Thus, if only the Coulomb-like potential was present, it would order 

the P-levels, in ascending order ^Pq, ^Pi, ^Pq. On the other hand, the existence solely of 

the scalar part would order them oppositely. The fact that the previous order is obtained 

suggests that the spin-orbit part due to the vector potential dominates, 

(iii). the ratios p (xfej(l^)) , P (Xw(l^)) i= (2++ - l++)/(l++ - 0++))can be calculated 

from Table ^. They are equal to 

piXbAlP)) =0.69 ±0.03 
p(Xw(2P)) =0.56 ±0.09 



which are in a very good agreement with experiment [41 



p{Xbj{lP)) =0.67 ±0.03 
p{Xbj{2P)) =0.58 ±0.03 

(iv) . the states are characterised by small speed (((f /c)^) < 0.11) which make the applica- 
tion of the Pauli approximation and perturbation theory possible, (^{v/c)'^'^ is proportional 
to the A of the corresponding state and inversely proportional to the mass of the bottom 
quark, 

(v) . ro depends on the energy of the state and tq ~ 12 fm. [ -Rinitiai, ''"o ) that is clearly 
larger than the range within which the interaction takes place. If we continue the integra- 
tion of the differential equations for bigger interparticle distances, we notice that at r = ro 
both the wavefunction and its derivative vanish due to confinement. 

(vi) . the splittings of the IS', 25', 35', 45' states are given from Table |^. 

The successful description of bottomium by means of the Breit equation can serve as a 
first indication that the equation can be used in short-range interactions, however, before 
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we generalise, it would be wise to study the spectrum of other quarkonia composed of 
lighter quarks. 



5.2 The Charmonium system 



In this section, a similar procedure is going to be followed in order to calculate the 
charmonium spectrum. The twelve states we will consider have the quantum numbers 
J^^ = 1 , 0^^, 2++, 0~+, 1^^. As in the previous subsection, the widths of the decays 
of the S states into e~^e~ will be employed and the relations between the As of the radial 
part of the wavefunctions will be eventually revealed. The decays lead to ratios given by 
the following expressions 



r{2S/lS) 
r{3S/lS) 
r{AS/lS) 





2 

V'V(25)(0) 




^J/V(1S)(0) 


^J/V(iS) 


2 




^j/V(is)(0) 


^J/V(iS) 


2 

^^(45) (0) 


-^?/.(4S) 


?/'j/V(is)(0) 



0.41 ±0.05 



0.14 ±0.03 



0.09 ±0.02 



(53) 
(54) 
(55) 



The parameters that need to be evaluated are three, the mass of the charm quark, = m, 
Xis and Xip. It is assumed that the coupling constant and string constant k are the 



same as in the bottomium system. This lies in the so-called flavour independence of 
the strong interactions which suggests that the interaction between c, c does not differ from 
that between b, b. The same holds true for any combination of two quarks. Since all quarks 
exist in the same three colour states, they must have identical strong interactions. 

The "fitting procedure" is going to be followed, in order to determine the three pa- 
rameters. We first impose the constraint that the mass of the charm quarks lies within the 
range 1 GeV <m< 1.6 GeV (running mass) The energies of the states J /ip{lS),ip{2S), 
Xci(l-P) are chosen as input values and the ratios [ (|53D -(^) as well as the differential equa- 
tions Eqs.[(|3ll)-(|33D] are used. The latter are integrated in the range [10~^°rc,3 fm] where 
Tc = a^/A. By solving the differential equations and by taking into account the relations 
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between the parameters imposed by [(p3D-(^)], we obtain the following results []T0 



mc= = 1.572 ± 0.009 GeV 

\is ~ 0.395 (GeV)2 (56) 
Aip - 0.074 (GeV)2 

which, in turn, give 

A25 0.083 (GeV)2 

A35 § 0.040 (GeV)2 (57) 
A45 ^ 0.030 (GeV)2 

Now that the parameters have been estimated, they can be used to calculate the rest of 
the energy levels of the system. The application of Eqs.[(^)-(^)] and Table |^ lead to the 
results summarised in Table ^. These values can be used to calculate the S and P splittings 
which have been measured in the laboratory and their discrepancy and, eventually, the 
success of the theory can be evaluated (Table ^ Table |10D . We notice there is a discrepancy 
of 18% between the theoretical and the experimental splitting of the Xc2(l-P) and Xci(l-P) 
states. The ratio p = (2++-1++)/ (1++ - 0++) equals p = 0.420+0.005 which is well below 
the experimentally calculated (= 0.49) Thus, although, some of the energy levels seem 
to be in a very good agreement with experiment, the results should be taken with great 
care. The reason why the charmonium system is not described entirely adequately by the 
theoretical model is that our theory considers the static interaction as the main contributor 
to the system's energy and it treats the retarded terms perturbatively. The problem is that 
for charmonium, this scheme does not provide exactly the correct energy as it nearly does 
for bottomium, resulting to larger deviations from experimental values. In contrast to 
bottomium where the corrections due to the Breit terms constitute a small proportion of 
the total energy, in charmonium they are quite large. Anyway, larger corrections to the 
static energy (which, also, contains relativistic contributions) were expected since these 
terms are inversely proportional to the mass of the constituent particles. By looking at the 
expectation value of the square of the speed of the systems (Table |^, Table ^ and Table 
H), we can conclude that, indeed, although the relativistic corrections are not negligible, to 
some extent, both bottomium and charmonium can be regarded as non-relativistic systems, 
with this assertion suiting more to bottomium. 
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Another reason for the failure to obtain the correct results is the form of the potential 
that is supposed to describe strong interactions. Although the potential employed exhibits 
the correct behaviour at small and large distances, we are not able to determine what 
it looks like in between. In other words, we do not know its exact form. Certainly, a 
knowledge of this would lead to better results. 

In addition, we should not forget the order to which the energy values has been calcu- 
lated. If the wavefunction used, corresponded to the solutions of the Coulomb problem, 
the spectrum would contain corrections of order up to a^. For the wavefunctions of the 
harmonic oscillator, the energies contain terms of order up to a^. If we wish to calculate 
the energies more accurately, we need to consider higher order terms since their inclusion 
would provide an essential contribution to the energy levels, however, again, the results 
should not agree exactly with experiment. This is because it is not possible to take into 
account all diagrams. Thus, we should better compare our results with other theories. 
We consider the study conducted by Eichten et al. (they, also, used a funnel potential 
as a possible form of the interaction) |]28| on charmonium. They fitted the parameters 
{us = 0.39, K, = 926 GeV/fm) and calculated the energies summarised in Table 0, Table 
P^ . The asterisk (*) refers to the states used to fit the parameters. Their values are slightly 
different from those we used because at the time those calculations were carried out, the 
experimental data differed from the recent results. The P- states, also, in these tables 
represent the centre-of-mass state of the j = 0, j = 1, j = 2 states. 

The coupling constant of the interaction and the QCD scale parameter were calculated 
in the previous subsection where we studied the bottomium spectrum. One might wonder 
whether it would be wise to allow the values of these parameters to be derived from the 
study of charmonium, instead. The problem is that in order to estimate these values, one 
approximation has been made already by optimising the radial wavefunctions of the three- 
dimensional harmonic oscillator (in addition to the approximation concerning the form of 
the potential and the approximation related to the effective interaction). A more rigorous 
and accurate procedure should assume the use of the wavefunctions of the funnel potential. 
It is true that charmonium exhibits a more relativistic behaviour than bottomium and 
a very careful consideration should be taken, otherwise there is a danger of obtaining 
unreliable results. The fact that bottomium consists of heavier particles and thereby the 
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system is regarded as less relativistic implies that such approximations as those performed 
throughout this work are less likely to result in unacceptable and meaningless results if the 
parameters are calculated for bottomium. 

6 Conclusions 

In this paper, the application of the Breit equation to bound state of systems charac- 
terised by short-range interactions was studied. Despite the clear lack of covariance of the 
equation, we assumed that it can be considered to be approximately Lorentz invariant for 
systems interacting through a very strong short-range potentials such as quarkonia and 
more specifically, bottomium and charmonium. A funnel potential was introduced as a 
candidate potential describing QCD and the Pauli approximation of the Breit equation 
was employed in order to distinguish the contribution of the static terms of the interaction 
from that part of the interaction responsible for retardation. The solution of the Breit 
equation helped to obtain the spectrum of bottomium and charmonium and to compare 
them with experimental data. The S and P states as well as their splitting were calculated 
and it turned out that the results are, in general, in a very good agreement with experi- 
ment. The Breit equation was, also, used to determine the coupling constant and the scale 
parameter of the strong interaction and they were estimated to be as = 0.36 ± 0.04 and 
A = 409 ± 40GeV, respectively 

A part of the success was attributed to the fact that the small speed of the states 
allowed the Breit terms to be treated perturbatively despite the relatively large coupling 
constant (a^/m^ is small while as is large). The large mass of the constituent particles 
makes the systems less relativistic, therefore the static part of the interaction contributes 
much more to the binding energy than the Breit terms. 

Unfortunately, for systems that are characterised by larger speed or, equivalently, by a 
smaller mass, the discrepancy between experiment and the theory is large enough to raise 
concerns over the reliability of the Breit equation as a satisfactory theory. This suggests 
that, in that case, the Breit equation is not adequate to describe systems of two fermions 
and an alternative equation should be employed. 
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Meson 


Energy due to 


Energy corrections due to 


States 


the Static potential 


/ — \ 
the Breit terms ^index = 


l{lb) 




1 Q n A 


1 [Z^) 


-ecstatic I J- [Zo)) 


— oU.U/i25 


T(3S) 


£'static(T(35')) 


-183.2^35 


T{AS) 


^static(T(45)) 


-279.7A4S 


T{5S) 


^static(T(55)) 


-371.0^55 


T{6S) 


i?static(T(65)) 


-458.1^5 






-320.0Au' 




£'static(X6l(lP)) 


-I6O.OA1P 


X62(1P) 


^static (Xfe2( IP)) 


-70.4Aip 


Xw{2P) 


Pstatic(X60(2P)) 


-5I8AA2P 


Xbii2P) 


Est.tUXbi{2P)) 


-3IOAA2P 


X62(2P) 


E,t.tic{Xb2{2P)) 


-I93.9A2P 



Table 1: Energy levels of bottomium as a function of the parameters ag, rub, As. 



Meson 




Energy due to 


Energy corrections due to 


States 




the Static potential 


the Breit terms |^Andex = ^^\/%^^ 




0-+ 


^static (^?( IS")) 


-48.0^15 


V{2S) 


0-+ 


£^static(^7(25)) 


-416.0^25 


r}{3S) 




^static (35")) 


-6O3.2A35 


77(45) 


0-+ 




-769.7^45 


h{lP) 


1+- 


£^static(/i(lP)) 


-128.0^iP 


h{2P) 


1+- 


E,t.tic{h{2P)) 


-268.8A2P 



Table 2: Energy levels of not experimentally measured bottomimum states as a function 
of the parameters cts, mb, As. 
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Meson 
otates 


jPC 


Energy 

(Experim.) 

(MeV) 


Energy (Theory) (MeV) 


1 

\u) / 


Static part 
(MeV) 


Retarded part 
(MeV) 


T{1S) 
2nd 


\— 


9460.37 ±0.21 


9460.37 ±0.21 * 


0.050 


9454.89 ± 0.21 


+5.48 ±0.01 


T{2S) 
8th 


\— 


10023.30 ± 0.31 


10023.30 ±0.31 * 


0.083 


10027.40 ±0.31 


-4.101 ±0.007 


T{3S) 
14th 


1~ 


10355.3 ±0.5 


10355.3 ±0.5 * 


0.111 


10362.9 ±0.5 


-7.42 ± 0.01 


T{4S) 
16th 


1-- 


10580.0 ±3.5 


10580.0 ±3.5 * 


0.087 


10584.9 ± 3.5 


-4.929 ±0.008 


T{5S) 
17th 


1~ 


10865 ± 8 


10882 ± 8 


0.125 


10890 ± 8 


-7.86 ±0.01 


T{6S) 
18th 


\— 


11019 ±8 


11037 ±9 


0.083 


11041 ± 9 


-3.944 ±0.007 


V{1S) 
1st 


0-+ 




9393.9 ±0.3 


0.050 


9426.8 ±0.3 


-32.88 ±0.08 


7j{2S) 
7th 






9973.3 ±0.4 


0.083 


9994.6 ± 0.4 


-21.32 ±0.04 


rj(3S) 
13th 


0-+ 




10324 ± 2 


0.111 


10348 ± 2 


-24.42 ±0.04 


15th 


0-+ 




10553 ± 1 


0.087 


10567 ± 1 


-13.56 ±0.02 


XwilP) 
3rd 


0++ 


9859.8 ±1.3 


9858.9 ±0.4 


0.043 


9869.1 ± 0.1 


-10.15 ±0.02 


Xi,i{lP) 
5th 


1++ 


9891.9 ±0.7 


9891.9 ±0.7 * 


0.043 


9897.0 ±0.7 


-5.077 ±0.001 


xUiP) 

6th 


2++ 


9913.2 ±0.6 


9914.7 ±0.9 


0.043 


9916.9 ±0.9 


-2.234 ± 0.004 


HIP) 
4th 


1+- 




9906 ± 1 


0.043 


9910 ± 1 


-4.062 ± 0.007 



Table 3: The energy spectrum of the bottomium system derived from experiment (third 
column), the Breit equation (fourth and fifth columns) and the expectation values of the 
square of the velocity operator v for each state (sixth column). The asterisk (*) refers to 
the states used to fit the parameters. The position of every state in the spectrum (for the 
states we have considered) is mentioned in the first column. 



34 



Meson 
States 




Encrgy 
(Experim.) 

(McV) 


Energy (Theory) (MeV) 




Static part 
(MeV) 


Retarded part 
(MeV) 


Xbo{2P) 
9th 


0++ 


10232.1 ±0.6 


10234.3 ±0.2 


0.087 


10253.9 ±0.2 


-19.60 ±0.03 


X6i(2P) 
11th 


1++ 


10255.2 ±0.5 


10255.2 ±0.5 * 


0.087 


10266.9 ±0.5 


-11.74 ±0.02 


Xb2{2P) 
12th 


2++ 


10268.5 ± 0.4 


10266.9 ±0.7 


0.087 


10274.3 ±0.7 


-7.33 ±0.01 


h{2P) 
10th 


1+- 




10260.2 ±0.9 


0.087 


10270.4 ± 0.9 


-10.17 ±0.02 



Table 4: The energy spectrum of the bottomium system derived from experiment (third 
column), the Breit equation (fourth and fifth columns) and the expectation values of the 
square of the velocity operator v for each state (sixth column). The asterisk (*) refers to 
the states used to fit the parameters. The position of every state in the spectrum (for the 
states we have considered) is mentioned in the first column. 



AM 


Theory (MeV) 


Experiment (MeV) 


Deviation (%) 


m(x6i(lP))-m(xbo(lP)) 


33.0 ±0.8 


32 ±1 


3.1 


m{xm{lP))-m{xbi{lP)) 


22.8 ± 1.1 


21.3 ±0.9 


7.0 


m(x6i(2P))-m(x6o(2P)) 


20.9 ±0.8 


23.1 ±0.8 


9.5 


m(x62(2P))-m(xi,i(2P)) 


11.7 ±0.9 


13.3 ±0.6 


12.0 



Table 5: Splittings of the IP and 2P levels of bottomium. 
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AM 


Theory (MeV) 


m{T{lS))-m{r]{lS)) 


66.5 ± 0.4 


m{T{2S))-m{r]{2S)) 


50.0 ±0.5 


m{T{3S))-m{r]{3S)) 


31 ±2 


m{T{AS))-m{r]{AS)) 


27 ±4 



Table 6: Splittings of the S levels of bottomium. 



Meson 


Energy due to 


Energy corrections due to 


States 


the Static potential 


the Breit terms ^^index ^ 


J /lb (IS) 


Estut 


ciJ/iHis)) 


+8.0Ai,9 


c{2S) 


-E'slal. 




-80.0-425 




-E'stat 


c(^(3S)) 


-183.2^35 




-^stat 


cms)) 


-279.7v445 


VcilS) 


E&tat 


MIS)) 


-48.0^15 


Vci2S) 


EstaX. 


M^s)) 


-416.0^25 




-E'stat 




-603.2^35 


Vci^S) 


-E'stat 


Mm 


-769.7^45 


hc{lP) 


Eiit-At 


MIS)) 


-128.0Ais 




Estat 


c(Xco(lP)) 


-32O.OA1P 


Xci(lP) 


-^stat 


c(Xcl(lP)) 


-I6O.OA1P 


Xc2(lP) 


Estat 


c(Xc2(lP)) 


-70.4Aip 


Xcl(lP) 


-E'stat 


c(Xcl(lP)) 


-I6O.OA1P 


Xc2(lP) 




c(Xc2(lP)) 


-70.4AiP 



Table 7: Energy levels of charmonium as a function of the parameters cts, mc, 
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Meson 
btates 




Energy 

(bxperim.j 

(MeV) 


Energy (Theory) (MeV) 




Static part 
(MeV) 


Retarded part 
(MeV) 


J/i^ilS) 
2nd 


1~ 


3096.88 ± 0.04 


3096.88 ±0.04 * 


0.24 


3078.68 ± 0.04 


+18.2 ±0.1 


i,{2S) 
Stli 


1— 


3686.00 ±0.09 


3686.00 ± 0.09 * 


0.47 


3703.58 ± 0.08 


-17.58 ±0.03 


^(35) 
10th 


1— 


4040 ± 10 


4106 ± 10 


0.36 


4120 ± 10 


-13.54 ±0.02 


^{AS) 
12th 


1— 


4415 ± 6 


4454 ± 8 


0.37 


4467 ± 8 


-13.44 ±0.02 


r^{lS) 
1st 


0-+ 


2979.8 ±2.1 


2987 ± 2 


0.24 


3095 ± 2 


-108.0 ±0.8 


r^{2S) 
7th 


0"+ 


3594 ± 5 


3601 ± 1 


0.47 


3692 ± 1 


-91.4 ±0.2 


7^(35) 

9th 


0-+ 




4023 ± 1 


0.36 


4068 ± 1 


-44.59 ± 0.08 


nth 


0-+ 




4379.1 ± 0.9 


0.37 


4416.1 ±0.9 


-36.98 ±0.07 


Xco(lP) 
3rd 


0++ 


3417.3 ±2.8 


3422 ± 1 


0.30 


3481 ± 1 


-59.0 ±0.1 


Xci(lP) 
5th 


1++ 


3510.53 ±0.12 


3510.53 ±0.12 * 


0.30 


3510.02 ± 0.12 


-29. 19 ±0.05 


Xc2{lP) 

6th 


2++ 


3556.17 ±0.13 


3547.9 ± 0.2 


0.30 


3560.9 ± 0.2 


-12.98 ±0.02 


h{lP) 
4th 


1+- 


3526.14 ±0.24 


3483 ± 2 


0.30 


3507 ± 2 


-23.60 ± 0.04 



Table 8: The energy spectrum of the charmonium system derived from the experiment 
(third column), the Breit equation (fourth and fifth columns) and the expectation values 
of the square of the velocity operator v for each state (sixth column). The asterisk (*) 
refers to the states used to fit the parameters. The position of every state in the spectrum 
(for the states we have considered) is mentioned in the first column. 
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AM 


Theory (MeV) 


Experiment (MeV) 


Deviation (%) 


miij{lS))-m{vc{lS)) 


110 ±2 


117 ±2 


6.0 


m{ij{2S)) -m{ric{2S)) 


85 ± 1 


92 ±5 


7.6 


m{^{3S))-m{7],{3S)) 


83± 10 






m(^(45)) -m(r;e(4^)) 


75 ±8 







Table 9: Splittings of the S levels of charmonium. 



AM 


Theory (MeV) 


Experiment (MeV) 


Deviation (%) 


m{xci{lP))-m{xco{lP)) 


89 ±1 


96 ± 1 


4.3 


m(xc2(lP))-m(xci(lP)) 


37.4 ± 0.2 


45.6 ±0.2 


18.0 



Table 10: Sphttings of the IP states of charmonium. 
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Charmonium rric = 1.84 GeV 


State 


Energy (MeV) 




3095* 


^P{2S) 


3684* 




4110 


^(45) 


4460 




4790 


XcilP) 


3522* 



Table 11: The charmonium spectrum according to Eichten et al. 



Bottomium rrih = 5.17 GeV 


State 


Energy (MeV) 


T{1S) 


9460* 


T{2S) 


10050 


T{3S) 


10400 


T(4^) 


10670 


T(5S) 


10920 


T{QS) 


11140 




9960 


Xb{2P) 


10310 




10600 



Table 12: The bottomium spectrum according to Eichten et al. 
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